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Abstract. We present the analysis of 65 hours of high speed photometric observations of HD 12098 taken from 
State Observatory, Naini Tal and Gurushikhar Observatory, Mt. Abu on sixteen nights spanning from November 
1999 to November 2000. HD 12098 is the first rapidly oscillating Ap star discovered from the 'Naini Tal-Cape 
survey for northern hemisphere roAp stars'. It is the 32nd in the complete list. HD 12098 exhibits one predominant 
mode of oscillation at = 2.1738 mHz. The second-most significant frequency in our data is at 1/2 ~ 2.1641 mHz 
with a 1 cycle/day alias ambiguity. We argue that 1/2 is a rotational sidelobe of ui, rather than an independent 
pulsation mode. Evidence for the presence of two other frequencies at 2.1807 and 2.3056 mHz is also presented. 

Key words. Stars: chemically peculiar - Stars: oscillations - Stars: individual: HD 12098 
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1. Introduction 

Many cool Ap stars exhibit anomalously strong Sr, Cr 
and Eu lines coupled with high magnetic fields frequently 
as strong as several kilogauss. A few of these chemically 
peculiar stars exhibit broadband light oscillations in the 
4-16 minute period range with amplitudes of few milli- 
magnitude (mmag). This subset of Ap stars is termed as 
rapidly oscillating Ap (roAp) stars and the oscillations 
are attributed to low-order (Z < 3), high overtone {n ^ I) 
non-radial pulsation modes (see Kurtz 1990, Martinez & 
Kurtz 1995 for reviews). Long term variations in ampli- 
tude of pulsations are observed in many roAp stars. This 
modulation can be explained in the light of the oblique 
pulsator model. In this model the pulsation and magnetic 
axes are oblique with respect to the rotation axis of the 
star. The modulation in the amplitude is then caused 
by the orientation effect due to the periodic variation 
in the angle between the pulsation axis and the line of 
sight. This effect seen in many roAp stars is therefore ex- 
tremely useful in determining the rotation period of the 
star (Kurtz [19821 ). 

In order to detect new roAp stars in the northern 
hemisphere, the "Naini Tal- Cape Survey" wa s init iated 
in 1997 (Ashoka et al. ^OOOt Martinez et al. |200lD , be- 
cause 28 of the previously known 31 roAp stars are in 



the southern hemisphere. HD 12098 is the first roAp 
star discovered from this survey. The co-ordinates of this 
star (aaooo = 02^*00™. 7, (Jaooo = 58° 31'. 6 ) make it the 
northern-most k nown roAp star. HD 12098 is of spectral 
type FO (Olsen 198C) and has an apparent magnitude 
mv = 7.974. It has Stromgren color indices [b — y) = 
0.191, mi ^ 0.3 28, ci = 0.517, [3 = 2.796 (Hauck & 
Mermilliod |1998[). bnix = -0.122 and (Sci = -0.279 



Send offprint requests to: S. Seetha (seetha@isac.ernet.in) 
* JAP, Dept. of physics, IISc, Bangalore-560 012 



are calculated using the calibrations of Crawford (1975). 
Recently Wade et al. (2001) have detected a variable lon- 
gitudinal magnetic field of « 2kG in HD 12098. They 
also predict a rotational period of the order of one week. 
HD 12098 was discovered in 1999 from Naini Tal as a 
roAp star having a period of 7.68 minute with an ev- 
idence for multi-periodicity (Martinez et al. 2000) and 
amplitude modulation. We obtained about 65 hours of 
high speed photometric data as part of follow-up obser- 
vations out of which around 45 hours are from Mt. Abu 
during October 2000. Analysis of these data indicates the 
presence of multi-periodic oscillations in the star. Among 
these, one frequency is most likely due to the rotational 
modulation of the predominant frequency. This supports 
the oblique pulsator model for this star as is the case with 
several other roAp stars. Probable evidence of the other 
rotational side-lobe and the presence of other modes of 
pulsation are presented here. Details of observations and 
data analysis are given in the next section while results, 
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Table 1. Journal of high speed photometric observations 
of HD 12098. The second column is the BJD at the center 
of the run, third column is the total duration of observa- 
tions. Last column lists the amplitude corresponding to 
the frequency 2.1738 mHz for the individual nights (re- 
fer sec. 2.3). 



Date 


Time(BJD) 


Duration 


Amp 




(2451000-f) 


(hours) 


(mmag) 


21 Nov 1999" 


504.27369 


2.51 


0.689 


20 Dec 1999" 


533.15988 


2.00 


0.381 


21 Dec 1999 


534.18171 


1.24 


1.461 


22 Dec 1999" 


535.17755 


3.12 


1.412 


08 Oct 2000" 


826.30268 


1.10 


0.370 


11 Oct 2000" 


829.27440 


2.64 


1.199 


12 Oct 2000" 


830.25841 


1.58 


0.834 


15 Oct 2000" 


833.28734 


2.01 


0.673 


16 Oct 2000" 


834.94355 


1.42 


0.800 


18 Oct 2000* 


836.39091 


5.19 


1.054 


19 Oct 2000* 


837.33851 


7.75 


0.567 


20 Oct 2000* 


838.32158 


7.76 


0.534 


21 Oct 2000* 


839.31472 


7.06 


0.657 


22 Oct 2000* 


840.29415 


7.82 


1.904 


23 Oct 2000* 


841.31576 


8.85 


1.196 


07 Nov 2000" 


856.21624 


2.64 


1.469 


Total duration 


64.99 hrs 





"Observed at Naini Tal 
'Observed at Mt. Abu 

discussion and conclusions are given in the remaining part 
of the paper. 

2. Observations and data analysis 

The details of observations, data reductions and pulsation 
analysis are given in the following subsections. 

2.1. Data collection and Reduction 

The pulsations in the star HD 12098 were discovered on 
the night of 21/22 November 1999 during the roAp star 
survey observations at the 1.04 m Sampurnanand tele- 
scope of state observatory, Naini Tal, India. Follow up 
observations were carried out at Naini Tal and also at the 
1.2 m telescope at Gurushikhar, Mt. Abu, India. The log 
of observations is given in Table 0. 

On both the telescopes we used a three channel 
photometer (Ashoka et al. [200 1| ). During the survey at 
Naini Tal, we often use only the main channel of the pho- 
tometer, while the follow up observations at Mt. Abu are 
conducted using all the three channels, monitoring a field 
star and sky simultaneously. HD 12098 is a double star 
with a companion which is 2.6 magnitude fainter and at 



an angular separation of 6" (Olsen 1980). During obser- 
vations, both stars are kept well within the diaphragm 
of 30". The data are acquired as continuous 10 s inte- 
grations in th e Joh nson B-filter using the Quilt-9 software 
(Nather et al. |1990D running on a PC. Altogether, data for 



HD 12098 are obtained on sixteen nights. Out of them, the 
nights from Naini Tal are with durations of observations 
ranging from 1.1 to 3.1 hours while the six nights' data 
from Mt. Abu are of relatively longer duration ranging 
from 7.06 to 8.85 hours except on one night where the du- 
ration is 5.19 hour. The data obtained at both the places 
are corrected for dead time, sky background and earth's 
atmospheric extinction trend. 

The data sets are in time expressed as barycentric 
Julian day (BJD) versus fractional intensity with respect 
to the mean intensity of the star. Since the fractional in- 
tensity variations are very small these values are given as 
mmag. Figure ^ presents typical light curve of the star 
HD 12098 obtained on 21/22 Oct 2000 and 22/23 Oct 
2000. The presence of the 7.6 minute oscillation and the 
day to day variation in pulsation amplitude is clearly ev- 
ident in the light curves. 
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Fig. 1. A representative light curve of the star HD 12098 
obtained on nights BJD 2451839 and 2451840. Each point 
represents the average of four 10 s integrations in the 
Johnson B filter. 



2.2. Pulsation analysis 

The presence of periodicity in the data-sets is searched for 
using the Discrete Fourier Transform (DFT) technique. 
All the individual runs exhibit a period around 7.7 min- 
utes with varying amplitudes. The data sets collected at 
Mt. Abu, all of which are of considerable length also ex- 
hibit the presence of a frequency at 2.17 mHz in each 
of the runs, but with night to night variation in ampli- 
tudes ranging from 0.5 to 1.9 mmag. This clearly indicates 
the likelihood of the star being multi-periodic. In order 
to resolve closely spaced frequencies, all the runs from 
Mt. Abu are concatenated together and subjected to the 
DFT. Figure || shows the amplitude spectrum of this data 
set in the 0-5 mHz range. The dominant frequencies are 
around 2.2 mHz and there is no excess power well above 
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the noise level in any other region of the frequency spec- 
trum up to the Nyquist frequency (50 mHz). There appear 
to be four frequencies present around 2.2 mHz which are 
given in column 1 and 2 of Table |^. 



Table 2. The dominant frequencies and amplitudes ob- 
tained for the Mt. Abu data. First two columns are the 
results of the Fourier transforms (FT) and the last two 
columns are obtained by non-linear least square fit to the 
data. 



FT 


NLSQ 


Frequency 
(mHz) 


Amplitude 
(mmag) 


Frequency 
(mHz) 


Amplitude 
(mmag) 


2.17372 
2.16421 
2.18060 
2.30559 


0.941 
0.487 
0.263 
0.276 


2.17385±0.00006 
2.16414±0.00010 
2.18069±0.00016 
2.30559±0.00014 


0.918±0.036 
0.511±0.036 
0.276±0.036 
0.279±0.035 



The top panel of Fig. y shows an expanded view of the 
amplitude spectrum around 2.2 mHz. The highest peak in 
this spectrum is at vi = 2.1737 mHz with its 1 day alias 
pattern (the spectral window of the data set is shown in 
the inset). There is also an indication of a second period 
at a slightly lower frequency which is evident as a double 
peak in vi and its left side alias. In order to see whether 
it is an independent frequency or not, we prewhitened the 
main frequency vi from the time domain and performed 
an FT again. The resultant amplitude spectrum of the 
residual data is shown in the second panel of Fig. ^. This 
clearly establishes the presence of a second frequency. The 
peak with highest amplitude in this pattern corresponds 
to a frequency 1^2 = 2.1641 mHz. However, the amplitude 
of its alias at 2.1759 mHz is comparable. In FTs noise 
can alter the amplitudes of the peaks. Hence we cannot 
unambiguously identify which of the above two is the real 
frequency, although the existence of a second frequency 
is evident. Our identification of V2 = 2.1641 mHz is thus 
subject to a 1/day alias ambiguity. We again prewhitened 
this frequency V2 from the data set and the amplitude 
spectrum of the residual data is shown in the bottom panel 
of Fig. ^. There is an indication of two different frequency 
groups in this final spectrum, one around v-^ = 2.181 mHz 
and the other at 1/4 — 2.306 mHz. Both of these frequencies 
are only marginally above the noise level and are identified 
only by the typical window pattern they exhibit. It may be 
noted however that while the spacing between vi and 1^2 is 
about 9.5 /xHz, the spacing between the central frequency 
vi and is around 7 /iHz. If 2.1759 mHz is identified with 
1/2 then the spacing between i>i and 1^2 is as small as 2.1 
fiRz. 

The Mt. Abu data set is subjected to a non-linear least 
square fit in order to determine the errors on the frequen- 
cies and amplitudes of the four frequencies determined 
from the FT. The results of this analysis are given in col- 
umn 3 and 4 of Table |[ 



2.3. Search for long term period 

The pulsation amplitudes of HD 12098 exhibit consider- 
able variation from night to night. This is also indicative of 
long term periods. In order to search for periodic modula- 
tion of the pulsation amplitude on a time scale of days, we 
determine a nightly amplitude of the frequency vi by per- 
forming a linear least square (llsq) fit to each of the light 
curves acquired during 1999-2000. The amplitudes thus 
determined are listed in the fourth column of Table 0. 

The resulting data-set of time vs. amplitude (column 2 
vs. column 4 of Table |l|) was then subjected to a discrete 
Fourier transform, non-linear least squares (nlsq) fit and 
a periodogram analysis using the folding technique. All 
these methods yielded around five peaks, with the high- 
est corresponding to a period of 1.22 day (0.82 cycles per 
day). This period though not completely convincing pro- 
vided best fit (with a variance 17% lower than the next 
best period) for the llsq fitted amplitudes from Table [| as 
shown in Fig. ^. 

3. Results and Discussion 

The pulsation frequencies determined above are discussed 
below in terms of pulsation models. 

3.1. Frequency splitting and rotational period 

From the analysis of Mt. Abu data, we have determined 
the presence of at least two dominant frequencies {vi = 
2.1738 and 1^2 = 2.1641 mHz) in our data set along with 
two other frequencies {1^3 = 2.1807 and 1^4 — 2.3056 mHz). 
Several roAp stars are multi-periodic and exhibit many 
modes (viz., Kurtz 1990). Hence the first hypothesis to 
test is whether or not these are individual modes. The 
temperature of this star as estimated from both Stromgren 
indices and the spectral type, is 7820 ± 30 K (Moon & 
Dworetski 1985). From the evolutionary models the char- 
acteristic spacing (Ai/) between alternate T and 'n' modes 
for Ap stars has been given by Heller & Kawaler (1988) 
and Shibahashi & Saio (1985) as a function of tempera- 
ture. If the two dominant frequencies in HD 12098 whose 
temperature is 7820 K, have to be independent modes with 
a frequency spacing of only 9.5 fiHz or less, the theoret- 
ical models would then imply a highly evolved star with 
a luminosity in excess of 2L0 and mass > 2.5 Mq. None 
of the observational parameters show that HD 12098 is 
much more evolved than the other roAp stars. Moreover 
the typical spacing of alternate modes in other roAp stars 
ranges from 20 ^Hz to 80 ^uHz (Kurtz, 1990). We therefore 
conclude that the two dominant frequencies under discus- 
sion are unlikely to be independent modes. 

An alternative explanation may be that, it is the 'm' 
splitting due to rotation in / = 1 state (Ledoux 1951) 
which gives rise to (2/ -t- 1) modes spaced in frequencies as 

i^m - l^o = (1 - Cni)mQ, (1) 

where Cni is a constant and Q is the rotational frequency. 
For HD 12098, we have no apriori estimate of the rota- 
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Fig. 2. The Fourier transform of the Mt. Abu data set obtained during BJD 2451836 - 841. It is evident that the 
dominant frequencies located around 2.2 mHz are well above the noise level. 



tional period. If we use the theoretical values of C„; in 
the range of 0.01 to 0.001 as given by Shibahashi & Saio 
(1985) in Equ. (1), a 9.5 /iHz frequency spacing implies a 
rotational period of around 1.21 to 1.22 day. In case the 
two frequencies are caused due to 'm' splitting, then the 
weak 2.1807 mHz peak could be the third of the expected 
triplet. Whether this hypothesis is valid for HD 12098, 
requires an independent determination of its rotation pe- 
riod. Incidentally, such a measurement has often led to 
this theory being disproved in the case of other roAp 
stars (Kurtz 199C) favouring the oblique pulsator model 
instead. 

The pulsation characteristics of roAp stars are well 
understood in terms of the oblique pulsator model 
(Kurtz |1982|) and we apply it to HD 12098. The rota- 



tional period expected from the 9.5 /iHz frequency split- 
ting of the two dominant frequency peaks is 1.22 day. In 
this case the third component of a conjectured 1 = 1 fre- 
quency triplet can be expected at 2.183 mHz. It may be 
noted that the unequal frequency spacing of the two ob- 
served side-lobes with respect to the main peak may be 
due to the presence of magnetic field. 



Kurtz et al. (|1990|) and Shibahashi & Takata ( [1993D 
have described a generalized oblique pulsator model in 
which the effects of both the magnetic field and rota- 



tion are considered (see also Dziembowski & Goode ( |1985 , 
1986| )). The asymmetry in the amplitudes of the triplet 



can be used to estimate the effect of magnetic field vis-a- 
vis the rotational effect by the two parameters, 



Pi 



A, 



= tan i tan (3 



(2) 



where A-i, Afj, A^i are amplitudes of the triplet corre- 
sponding to frequency in increasing order, (3 is the angle 
between the rotation and magnetic axis (pulsation axis), i 
is the inclination angle of the rotation axis to the observer; 
and 



Pi 



A 



+1 



A-i 



CQ 



A 



+1 



(l)mag 



(l)mag 



(3) 



where C is a constant, fl is the rotation period and the 
perturbations due to the magnetic field depend on \m\ 

(l)mag 



such that uj — tj*^"^ + lu]^^^"^'^'^ . From the values given in 
Table |, we calculate Pi = 0.86 and P2 = -0.29. 

If vi and 1/2 are interpreted as components of a ro- 
tational multiplet in HD 12098, this would correspond 
to a rotation period of 1.22 day. This would then make 
HD 12098 the fastest rotating known roAp star. Assuming 
a radius of 2Rq the equatorial velocity is calculated to be 
w 80 km s^^. For Ap stars the limit of vs'mi is of the 



order of 10 km s"! (Wolff |1983| ). Hence for HD 12098, 
this would imply an i « 10° and f3 « 78°. Based on the 
magnetic field studies, observations of pole reversal would 
indicate j3 > 80°. However, we remind the reader that 1^2 is 
subject to a 1 cycle/day alias ambiguity. In the event that 
1^2 is identified with 2.1759 mHz instead of 2.1641 mHz, 
the spacing between i>i and 1^2 would be 2.1 fjRz. This 
would imply a rotational period of 5.51 day. Our present 
data cannot distinguish between these possibilities, and 
the issue will only be resolved by acquiring more data. 
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Fig. 3. A close-up view of the FT of the data set BJD 2451836-841. The second panel is the amplitude spectrum of the 
data after subtracting the frequency i>i — 2.1738 mHz from the time domain. The third panel is amplitude spectrum 
after subtracting the frequencies vi — 2.1738 and 1^2 — 2.1641 mHz from the time domain. The inset of panel I is the 
spectral window for the data-set calculated by sampling a noise-free sinusoidal with the frequency and amplitude of 
the highest peak in Fig. ||. 



3.2. Other modes 

The other frequency at the limit of detection in our com- 
plete data set is at 2.3056 mHz. This peak also exhibits 
the typical alias pattern and is comparable in amplitude 
to the 2.1807 mHz peak. It may be noted however that this 
frequency has been detected in some of our earlier data- 
sets. The runs in which this frequency is predominantly 
above the noise level corresponds to observations on BJD 
2451504.27 and 533.16. Hence this may be either a short 
lived mode or a mode which is severely modulated. In ad- 
dition we have also detected excess power around 4 mHz 
on several nights. 

3.3. HD 12098 and other roAp stars. 

HD 12098 appears to exhibit a simple power spectrum 
at first glance. There are several roAp stars such as 



HD 196470 and HD 193756 (Martinez et al. |l99lD which 
exhibit only a single frequency. However, HD 12098 ap- 
pears to have few modes exhibiting effects of rotation in 
the star. Since the amplitude of the main pulsation does 
not exceed 2 mmag, it is one of the low amplitude roAp 
pulsator. Further observations may yield other low ampli- 



tude modes in the star, if there are any. If the rotation 
period of the star is 1.22 day or close to it, then it be- 
comes the fastest rotating roAp star. HD 6532 whi ch has 
a rotational period of 1.9455 day (Kurtz & Marang 1987) 
also has a similar amplitude spectrum with six frequency 
peaks (Kurtz & Cropper 1987] ). Anot her r oAp star which 
is very similar is HD 80316 (Kurtz 1990). The rotation 
period of this star is not yet confirmed, but is probably 
around 2 days. 

The estimated Smi (= -0.122) for HD 12098 is at the 
extreme edge of the selection criteria which sets the limits 
for probable roAp stars for this index (Martinez & Kurtz 
1995). However there are other roAp stars with even more 
negative 6mi (e.g., HD 101065). In this respect it may be 
noted that HD 80316 has 6mi = -0.118. The principal 
oscillation frequency of this star is 2.2516 mHz with an 
amplitude of 0.44 mmag. Hence HD 80316 can be consid- 
ered very similar to HD 12098 in several observed prop- 
erties including the frequency spectrum and the pulsation 
amplitude. 
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Fig. 4. Plot of pulsation amplitude of HD 12098 vs. BJD. X-axis is the BJD at the middle of the run and Y-axis is the 

pulsational amplitude in milli-magnitudc. The points arc the experimental results with the errors in the amplitude. 
The dashed curve is the sinusoidal corresponding to the period 1.22 day. 



4. Conclusions 

We have presented the photometric observations of a 
new roAp star HD 12098 discovered in the "Naini Tal- 
Cape survey for northern hemisphere roAp stars" . The 
detailed analysis of the observations shows the presence 
of two predominant frequencies. The main frequency is 
at 2.1738 mHz corresponding to a period of 7.67 minutes 
with an amplitude of 0.918 mmag. The second frequency 
can be either 2.1641 or 2.1759 mHz. This is most likely 
the rotational side-lobe of the main frequency vi . 

Hence we propose that HD 12098 is an oblique pul- 
sator which pulsates in the I = 1 mode with a main pe- 
riod of 7.67 minute. The rotation period of this star is of 
the order of a few days, but with our data we are unable 
to discriminate between possible rotation periods of 1.2 d 
or 5.5 d. Independent determination of the rotational pe- 
riod using differential photometry, or additional multi-site 
high-speed observations are required to resolve this ambi- 
guity and to improve our understanding of the pulsation 
spectrum of HD 12098. 
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